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Abstract: Entanglement (Renyi) entropies of spatial regions are a useful tool for character- 
izing the ground states of quantum field theories. In this paper we investigate the extent to 
which these are universal quantities for a given theory, and to which they distinguish differ- 
ent theories, by comparing the entanglement spectra of the massless Dirac fermion and the 
compact free boson in two dimensions. We show that the calculation of Renyi entropies via 
the replica trick for any orbifold theory includes a sum over orbifold twists on all cycles. In a 
modular-invariant theory of fermions, this amounts to a sum over spin structures. The result 
is that the Renyi entropies respect the standard Bose-Fermi duality. Next, we investigate 
the entanglement spectrum for the Dirac fermion without a sum over spin structures, and for 
the compact boson at the self-dual radius. These are not equivalent theories; nonetheless, 
we find that (1) their second Renyi entropies agree for any number of intervals, (2) their 
full entanglement spectra agree for two intervals, and (3) the spectrum generically disagrees 
otherwise. These results follow from the equality of the partition functions of the two theories 
on any Riemann surface with imaginary period matrix. We also exhibit a map between the 
operators of the theories that preserves scaling dimensions (but not spins), as well as OPEs 
and correlators of operators placed on the real line. All of these coincidences can be traced 
to the fact that the momentum lattice for the bosonized fermion is related to that of the 
self-dual boson by a 45° rotation that mixes left- and right-movers. 
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1. Introduction 



The quantum entanglement between spatial regions, as quantified by entanglement entropies 
and entanglement Renyi entropies, is an important tool for characterizing the infrared be- 
havior of extended quantum systems. In theories with a mass gap, these quantities provide 
a characterization of topological phases where no local order parameter exists [1, 2, 3, 4]. 
When the infrared limit is a nontrivial two-dimensional conformal field theory, computing 
the entropy of an interval (for example in a lattice model) provides an efficient way to de- 
termine its central charge [5] . The entropies of more than one interval in a CFT depend on 
the full operator spectrum, and therefore give more refined information about the theory; 
as an example, for a free compact boson the Renyi entropies of two intervals depend on the 
compactification radius [6, 7]. 

An important question is thus whether the entanglement entropies (by which we mean 
both von Neumann and Renyi entropies) of spatial regions in the ground state of a field theory 
characterize that theory, in the following precise senses: 

1. The entanglement entropies should be the same regardless of the presentation of the 
theory (up to non- universal cutoff-dependent terms). That is, the same theory could 
have two different Lagrangian descriptions, but the underlying spectra of states and 
local operators are the same; in this case, for a quantity to be universal it should give 
the same answer for both presentations. For example, in the case of the compact boson, 
the Renyi entropies are T-duality invariant [6, 7]. 

2. The entanglement entropies should distinguish different theories. 

Surprisingly, these statements appear to be challenged already by some very simple quantum 
field theories, namely the free massless Dirac fermion and the compact free boson in two 
dimensions. In the case of the Dirac fermion, the Renyi entropies for any number of intervals 
have been computed by Casini, Fosco, and Huerta [8, 9, 10]. In the case of the compact boson, 
the Renyi entropies for two intervals have been computed by Calabrese, Cardy, and Tonni 
[6, 7]. Since the boson at radius R = \^2R s d (where R S( ± is the self-dual radius) is known to be 
dual to a theory of a Dirac fermion, one might expect that the entropies computed by CFH 
would agree with those computed by CCT at that value of R. In fact, they do not, seeming 
to violate point (1) above. Various explanations have been put forward for this discrepancy, 
such as that the bosonization relating the two theories is a non-local transformation on the 
fields appearing in the path integral, and hence might not preserve the factorization of the 
Hilbert space according to spatial regions [10, 7]. Even more curiously, the Dirac fermion 
Renyi entropies do agree with those for the boson at R = R s d, despite the fact that these two 
theories are certainly not dual to each other. By conformal invariance the Renyi entropies 
for two intervals are effectively functions of the cross-ratio of the four endpoints, so this is 
a non-trivial agreement between an infinite number of functions of one real variable. This 
coincidence would seem to threaten point (2). 
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In this paper we will show that in fact both points continue to hold. The key to resolving 
the first discrepancy is to recall that there are two versions of the free Dirac fermion theory. 
We will show in Section 2 that the calculation of CFH [8] , as well as a subsequent calculation 
by Casini and Huerta giving the same result by a different method [9], specifically produces 
the entanglement entropies for the Dirac fermion theory without any projection on fermion 
number (and containing only NS-NS sector operators). This is not a modular-invariant theory, 
and is not dual to the boson at radius \f2R s &. Instead, it is only after a certain Z2 gauging, 
which introduces R-R operators and removes the fermionic ones, that the theory becomes 
modular-invariant and dual to the \/2i? s d boson. We will show that when the Renyi entropies 
of theories with such discrete gaugings (including orbifold theories) are calculated in terms 
of the partition functions on Riemann surfaces, one must perform a sum over over all twists 
by the gauge symmetry. This guarantees that the Renyi entropies for arbitrary numbers of 
intervals are invariant under Bose-Fermi duality. Thus, point (1) above is satisfied in this 
case. 

Having dispensed with the boson at R = \/2R s d, in the remaining sections we will explore 
the relationship between the original (unprojected) Dirac fermion and the self-dual boson. 
Our goals will be to understand the origin of the surprising agreement between their Renyi 
entropies for two intervals, to discover whether it extends to more than two intervals, and to 
see whether such a coincidence could happen in other (perhaps more complicated) theories. 
In doing so we will find that the theories are related by a new kind of quasi-duality, which 
we call a "real duality", that goes well beyond Renyi entropies. 

In Section 3 we will study the partition functions of the theories on Riemann surfaces, 
finding that they agree precisely when its period matrix is imaginary. Using the symmetries 
of the Riemann surface involved in computing the nth Renyi entropy for N intervals, we show 
that this condition holds when N = 2 for any n (explaining the agreement found before), and 
also for any N when n = 2, but not more generally; the Renyis for n > 2, N > 2 do indeed 
distinguish between these two theories, so condition (2) above is satisfied. (Therefore the von 
Neumann entropies, which are related to the Renyi entropies by an analytic continuation in 
n, also presumably distinguish between the theories for N > 2.) Since the Renyi entropies 
are known for the Dirac fermion for all iV [10], as a bonus of our analysis we learn what the 
n = 2 Renyis are for all N (see equation (2.3)). 

The agreement between the partition functions for imaginary period matrices is due to the 
following relationship between the theories: If we bosonize the Dirac fermion, its momentum 
lattice (which is simply Z 2 ) is related to the one for the self-dual boson by a 45° rotation (see 
figure 4). This rotation preserves the scaling dimensions of the corresponding momentum 
operators, but, since it mixes left- and right-movers, it changes their spins. Since the two 
theories also have the same oscillator structure, they have the same total spectrum of scaling 
dimensions, and hence the same torus partition function for imaginary r, i.e. on a rectangular 
torus. The agreement for higher-genus Riemann surfaces with imaginary period matrices is 
a generalization of this statement. 

In Section 4, we will use the 45° rotation on the momentum lattices to define a canonical 
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one-to-one correspondence between the operators of the two theories that preserves not only 
the scaling dimensions, but also (1) the OPE of any two operators that are separated by a real 
interval; (2) the correlator on the plane of arbitrary operators with positions on the real axis; 
and (3) the action of the mixed Virasoro generators L n + L n . We refer to this relationship 
between the two theories as a "real duality" . The coincidence of correlators is directly related 
to the statement about Riemann surfaces with imaginary period matrices, assuming a certain 
conjecture about their Schottky parameters. 

In Section 5 we extend our results on the agreements and disagreements between the boson 
and fermion theories to finite temperature and finite volume, and we discuss generalizations to 
other pairs of theories. We then discuss various larger issues connected to our work, returning 
in particular to the question we started with, whether entanglement entropies characterize 
theories. 

There are also three supplemental Appendices. Appendix A is a derivation of a specific 
expression for the partition function of the self-dual boson with radius R = 1, filling in a 
calculation needed in section 3. Appendix B summarizes the irreducible representations of 
the dihedral group D n , also needed for section 3. Finally, section C presents a rederivation 
of the results of [7] for the Renyi entropies for two intervals of a boson at any radius, and 
extends this result to the second Renyi entropy for any number of intervals. 

2. Entanglement entropies and discrete gauge symmetries for 1 + 1 CFTs 

In this section we will revisit a puzzle that arises in the computation of the entanglement 
entropies for the free massless fermion and free massless boson in 1 + 1 dimensions. We will 
begin in subsection 2.1, with a review of the definition of entanglement entropies in conformal 
field theory. In subsection 2.2, we will review recent results for the free Dirac fermion and the 
boson in 1 + 1 dimensions, which seem to indicate that the Renyi entropies are not invariant 
under Bose-Fermi duality [10, 7]. To prepare for a deeper investigation of this question, 
we will review the "replica trick" calculation of entanglement entropies in subsection 2.3. 
In subsection 2.4 we will extend that prescription to theories with discrete abelian gauge 
symmetries, which includes both bosonic orbifolds and modular-invariant fermionic CFTs. 
The invariance of the Renyi entropies under bosonization will follow automatically from the 
results of this subsection and old results about partition functions and bosonization [11, 12]. 
Finally, in subsection 2.5 we will resolve the puzzle by showing that the fermionic theory in 
question is not in fact dual to the free boson in question. 

The results we will review in subsection 2.3 also include an equality between the entan- 
glement Renyi entropies for two theories that are not dual to each other. Section 3 of this 
paper will be devoted to explaining this puzzling equality. 

2.1 Definition of entanglement entropies 

We will consider 1 + 1-dimensional conformal field theories (CFTs) on the plane. For thorough 
reviews of the entanglement entropies in these theories, we recommend [10, 7]. 
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Let A be the union of N intervals on the real axis, A = Vjf =l Ai, Ai = [ui, Vj\ (u\ < v\ < 
U2 < • • • )• We denote the complement as B = \jf =0 Bi, where Bi = [vi, Ui+i], with vq = — oo, 
u n+i = oo. (See Fig. 1.) We will assume that for a local theory, we can decompose the Hilbert 
space as H = Ha ®T~Lb, where Ha,b is the Hilbert space of degrees of freedom localized on 
A, B. In practice one must regularize the theory; if we put the theory on a lattice, then this 
decomposition should make sense. 



B A , B , A 2 B 2 A 3 B 3 

u, v, u 2 v 2 u 3 v 3 



Figure 1: A line divided up into consecutive intervals. We will be considering the density matrix for 
the field theory on the intervals A:=i,2,3 upon tracing out the local degrees of freedom in the intervals 

Bi=0, 1,2,3- 

The density matrix for the vacuum is p = |0)(0|; the reduced density matrix on A is 
PA = tr% s p. Its Renyi entropy of order n, called the "entanglement Renyi entropy" since it 
is a measure of the amount of entanglement between Ha and He, is defined by 

S n (A) = -*—lntr p n A , (2.1) 

1—71 

where n is a positive real parameter not equal to 1. Typically the Renyi entropy is computed 
for integer values of n. Knowing S n (A) for all integer n > 1 is then enough in principle to 
fix, by analytic continuation, the value for all positive real n. In particular, the value of the 
analytically continued function at n = 1 is the entanglement (von Neumann) entropy S(A) = 
— tr pa In pa- One can also consider various interesting linear combinations, such as the 
mutual Renyi information between two intervals: I n (A\ : A2) = S n (Ai)+S n (A2)— S n (AiUA2) 
(or, more generally, between two disjoint sets of intervals). In any computation, the Renyi 
and von Neumann entropies will diverge. The goal is then to extract universal, regulator- 
independent quantities. For example, the divergent parts of the Renyi entropies cancel for 
the mutual Renyi information 

One way to compute S n (A) is to find an explicit expression for the reduced density matrix 
PA in some basis, and from it directly compute trp^. To our knowledge the only theory for 
which this has been accomplished is the free massless fermion [9, 10]. The more common 
method, which we will review in §2.3, is the so-called replica trick, in which tr p n A is expressed 
in terms of the Euclidean partition function on an n-sheeted Riemann surface with branch 
cuts along the intervals Ai. 

2.2 Results for free CFTs 

The classic result by Holzhey, Larsen, and Wilczek, derived using the replica trick, is that 
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the Renyi entropies for one interval are the same for all CFTs, up to an overall factor of the 
central charge [13]: 

S n (A 1 ) = l(l+ 1 -)ln V -^, (2.2) 



6 \ n / 
where e is an ultraviolet cutoff length. 

On the other hand, for more than one interval the Renyi entropies depend on more than 
just the central charge. The only theory for which the entropies have been computed exactly 
for any number of intervals is the free Dirac fermion, which was accomplished both using the 
replica trick (by Casini, Fosco, and Huerta [8]) and by deriving an explicit formula for pa (by 
Casini and Huerta [9]). The result is a remarkably simple formula, in which the n-dependence 
factors out entirely from the dependence on the configuration of intervals: 1 

S^(A) = Ul + ±) E(A), E(A) = ^ln|t; J -n i |-^ln(n i -^)-^ln(^-^)-A r lne. 

i,j i<j i<j 

(2.3) 

Calabrese, Cardy, and Tonni [6] computed the Renyi entropies for two intervals for the 
compact boson at arbitrary radius R using the replica trick. Their result is quite a bit more 
complicated than (2.3), but it can be conveniently written as (2.3) plus a correction term: 

S*{A) = S®{A) + ^— ln.F*(s) • (2.4) 
1 — n 

The correction term is finite (does not involve e), and depends only on the conformally 
invariant cross-ratio of the four endpoints: 

(u 2 - ui)(v 2 - Vi) 
J-jf is a ratio of Riemann-Siegel theta functions 

(2.6) 



where rj = R 2 /R^ d (R s d is the self-dual radius), T is an x-dependent (n — 1) x (n — 1) matrix 
with the following entries: 

2%^. fTrk\ f2irk(r-s)\ 2 F 1 (k/n,l-k/n;l,l-x) 

1„ = — > sin — cos - . — r — , 2.7 

\n J V n J 2Fi(k/n,l-k/n;l,x) v 1 

and the Riemann-Siegel theta function (at the origin) is 

tf(0|r) = Yl e inTrsmrmS . (2.8) 

Although the expression for S^(A) is complicated, three key points are clear just from 

(2-6): 

lr The function 3 can also be expressed as S,(A) = — In | det M\, where M is an iV x TV matrix with entries 
Mij = e/(vj — Ui). The result (2.3) has a number of interesting properties; for example, the mutual Renyi 
insformations obey I„(A : B U C) = I n (A : B) + I„(A : C) for any three sets A, B, C of intervals [8]. 
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• Fn{x), and hence S^(A), are invariant under T-duality. So, at least in this case, the 
presentation of the theory does not affect the entanglement spectrum. 

• At the self-dual radius (77 = 1), J"^ sd (x) = 1 identically, hence S^ sd (A) = Sip (A) for 
any n and any A. 

• At R = V2R s d (v = 2 ), where the boson is dual to a theory of a Dirac fermion, 
J~n( x ) 7^ 1 (this can easily be confirmed numerically, in case the reader is worried 
about theta- function conspiracies), hence S^ Rsd (A) ^ Sn\A). 

The discrepancy between gV2Rsd(A) and s£\A) has been noted in the literature but has 
not been satisfactorily resolved. Casini and Huerta propose that the mismatch is due to the 
fact that the bosonization transformation is non-local; therefore, although the two theories 
have the same Hilbert space, the way that that Hilbert space gets cut up according to spatial 
regions in the two presentations might be different [10]. Calabrese and Cardy imply that the 
discrepancy is related to the Lagrangian used in computing the entanglement entropies of the 
fermion [7]. This is not an unreasonable thing to expect — at the level of the path integral, 
bosonization is not a local transformation of the fields we integrate over. It is therefore fair 
to ask whether it is a local transformation at the level of the Hilbert space, that is, whether 
the factorization of the Hilbert space by spatial regions is invariant under arbitrary duality 
transformations. 

However, interpreting the mismatch requires some care. The correct Bose-Fermi equiva- 
lence is between the boson at radius y/2R s d and the Dirac fermion gauged in a specific way by 
the Z2 fermion number [14, 11]. For example, the fermionic theory without such a gauging is 
not a modular-invariant theory, while the free boson is. 2 After this gauging, the spectra and 
the algebras of local operators are identical. Since this data defines a two-dimensional CFT, 
we might expect that computations of the position-space entanglement entropies should be 
the same whether computed in the bosonic or fermionic representation. More precisely, the 
scheme used for cutting off the theory may depend on the representation of the theory, but 
universal quantities such as the mutual information I n (A\ : A2) should not. 

In the remainder of this section we will explain this apparent mismatch. 

2.3 Renyi entropies via the replica trick 

Let us review the calculation of the density matrix p and the Renyi entropies via path integrals 
in a two-dimensional conformal field theory C (again, see also [7]). denotes all microscopic 
fields in the theory; a matrix element of p is p(4>i,(f>2) = (</>i|0)(0|(fo)j where | ^1,2) are field 
eigenstates with eigenfunctions (j>i^(x). These inner products can be represented via path 

2 Note that one may add additional "topological terms" to the bosonic theory which spoil modular invariance 
and lead to a theory which is precisely equivalent to a fermionic theory with fixed spin structure [12]. This 
cannot be the theory of the modular-invariant bosonic at R = i? a( i. 
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integrals, 

Pc/>exp - / dr da £(» . (2.9) 

•J — oo <J — oo 

0(x,O)=(/)i(x) 

To trace out the spatial region B, let us use 4>a,b to denote the function 4>{x 6 (yl,i3),r). 
We decompose the boundary conditions in (2.9) into those at x G A and x £ B, so that 

rO 



P(/>exp 



t/>(x£A,0)=<l> A (x) 
*(i6B,0)=^ B (i) 



dr / da C{<f>) 



(2.10) 



Tracing out the region B in the unorbifolded theory is now simple: for functions ^ 2 on A, 



Pa(<I>i, fa) 



1 



|0)(0|^ 



^2,<PB/ 



(2.11) 



where is the partition function of the CFT on the Riemann sphere. (See figure 2A.) The 
Renyi entropies are: 



1 f N 

1 ^ a=l 



(2.12) 



(A) (B) 




Figure 2: A: The reduced density matrix p(</>i, </>2) computed via a path integral on the complex 
plane with cuts Ai on the real line, boundary conditions <fi = (f>i at the top of the cut and cf> — 4>2 on 
the bottom of the cut. B: The Riemann surface S„jy constructed as an n-fold branched cover of the 
complex plane, with branch cuts at Ai glued together as shown. 

The integral in (2.12) can be done by "replicating" the Euclidean spacetime. The &th 
integrand in the product, pa(4>o,, 4>a+i), is the path integral on the complex plane, with cuts 
on the real line at Ai, and boundary conditions 4> a at the "bottom" of the cut and </> a +i at 
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the "top" of the cut. In taking the product and integrating over all of the a s, we are taking 
n copies of the plane and gluing them together in cyclic order by identifying the top of the 
cut on sheet k = 1, . . . n — 1 with the bottom of the cut on sheet k + 1, and glued the top of 
the cuts on sheet n with the bottom of the cuts on sheet 1. (See figure 2B). This is a singular 
Riemann surface £ ni jv with genus g = (n — 1)(N — 1), described as an n-fold branched cover 
of the sphere over N branch cuts. The result is that 

trpl=%^ (2.13) 

where Z n ^ is the partition function of the CFT on S nj jy- 

While we used a basis of field eigenstates to construct the theory, this is not necessary in 
principle. We could have used any other basis that respects the decomposition % = Ha^T^b- 

2.4 Including discrete gauge symmetries 

Next, consider orbifolds of our theory C by a discrete symmetry group G. For the sake of 
simplicity, let us consider the case of a bosonic orbifold C/G with the unorbifolded CFT C 
described as above by scalar fields </>, and a G-action 4> l— ^ 94 1 - The G-action could be a finite 
rotation or a discrete translation. We will restrict to the case that G is abelian. Before tracing 
out any spatial regions, (2.10) still holds. 

Breaking up 4>b into its values {4>Bi} on each interval Bi, the reduced density matrix for 
the intervals Ai in the orbifold theory is: 

PA(<k,<h) = Y Yl f^B 1 ---V<j )BN ( ( j> A = <p 1 ,{g i( j )Bi }\0)(0\ ( l> A = ( l> 2 ,{g i( l> Bi }). (2.14) 
1 {<?i}eG* 

Here gi £ G; in taking the trace over degrees of freedom in the intervals Bi, we have identified 
<j> up to discrete gauge transformations, so that the trace is being taken in C/G. Z\ is again 
just the path integral on the Riemann sphere, and its presence ensures that XxaPa = 1- The 
result is that the reduced density matrix is the sum over path integrals on the cut plane shown 
in figure 3, with each element of the sum corresponding to twists of the field by hi = g%gi-\ € G 
as one transports the fields around the cuts Ai. 

The fields on either side of all of the cuts are still untwisted as one transports them from 
— oo to oo along a curve parallel to the real line: they can be deformed along the imaginary 
axis to ±ioo, and the system is in the vacuum which is generally an untwisted state. Thus, if 
we take the sum of a left-directed contour above the cuts and a right-directed contour below 
the cuts, we can deform them so that they become the sum of single contours around each 
cut. Since the fields have zero twist around the initial contours, the products of the twists 
around all of the cuts must be equal to the identity: Y\a=i hi = 1. 

Next, we wish to compute the Renyi entropies, by computing 

t *PN= E fWd^A^pig^Ai),^^)) (2.15) 
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Once again we have sewed together fields fa along intervals A\ up to identification by the 
gauge group G. The result is 

e d-)«.W = tr^ = ^. (2.16) 

The expressions (2. 14), (2. 15) make it clear that Z n ^ is the partition function for the orbifold 
CFT on S nj jv, in which we have summed over all G twists about all non-contractible cycles 
(again, the cycles at infinity on each sheet are taken to be trivial). In other words, one treats 
the Riemann surfaces exactly as one would Riemann surfaces for string perturbation theory 
for orbifold backgrounds. 

We discussed bosonic orbifolds for simplicity, 
but this argument will be identical for gaugings 
of fermionic theories. When the fermion number 
itself is gauged, Z n ^ will correspond to the sum 
over all spin structures of fermion partition func- 
tions on £jv !n . 

We are now in a position to argue that the 
entanglement entropies are invariant under Bose- 
Fermi duality. This duality is between the mass- 
less Dirac fermion gauged by Z2 fermion number 
and the free boson on a target space circle with 
radius R = ^/2R S &. It is known that the equiv- 
alence holds for partition functions on any Rie- 
mann surface, so long as one correctly sums over 
the fermion spin structures [14, 11, 12]. Since the 
Renyi entropies are determined by these partition 
functions, they are guaranteed to match. 

More generally, the results of this section al- 
low one to discuss entanglement entropies for a large class of orbifold theories. There are two 
more complicated generalizations which we leave for future work. One is the case of orbifolds 
with discrete torsion, in which the different twisted partition functions on a Riemann surface 
are added with nontrivial phases. The second is the case of nonabelian orbifold groups. We 
suspect that the result will be the same — one treats the Riemann surfaces in the Renyi en- 
tropy calculations precisely as one would the Riemann surfaces in string perturbation theory 
calculations (but without the integration over moduli). 

2.5 Boundary conditions in the Dirac fermion calculations 

Given the results of the prior section, a candidate explanation for the discrepancy described 
in [10] is that they are working with a different gauging or with the ungauged theory. We 
will argue that the latter is in fact the case, by examining both the replica calculation [8, 10] 
and the direct construction of the reduced density matrix [9, 10]. 




Figure 3: The cut plane used to calculate 
the density matrix for 2 intervals. The fields 
on each side of the slit are independent, as 
in (2.14). For Renyi entropies of the reduced 
density matrix corresponding to the CFT 
vacuum, the fields are untwisted for a cir- 
cle which encloses all cuts A4. These can be 
deformed to the sum of the two loops shown 
which encircle the cuts. The bosonic fields 
can be twisted about these loops, so long as 
the ordered product is the identity. 
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We begin with the calculation of Z n ^ using the replica trick [8, 10]. In applying the 
replica trick to a theory containing fermionic fields, one needs to be careful about boundary 
conditions for the fermions. In a theory in which fermion number is gauged, it follows from 
the discussion in the previous subsection that one should sum over the partition functions 
with NS and R boundary conditions around all cycles. However, in an ungauged theory, 
there is a specific set of boundary conditions implied by the replica trick (just as the thermal 
partition function is computed with a specific boundary condition — namely antiperiodic — 
around the Euclidean time circle). We will not review the derivation here; it is given in [8] 
(below equation (6); see footnote 6 of [15] for an alternate derivation), but the result is that 
for even n one must include a sign-flip along the cuts that connect the first sheet to the 
last one. (These boundary conditions imply that, when passing to a single-valued coordinate 
system in the neighborhood of a branch cut, there is no spin field inserted at the branch 
point. They also imply NS boundary conditions on all the basis cycles we will use in Section 
3, shown in figure 5.) Since no sum over boundary conditions was performed, it is clear that 
the calculation in [8] is done in the ungauged theory. Indeed, in the case N = n = 2, where the 
replicated surface is a torus, one can directly reproduce the CFH result from the well-known 
torus partition function for a fermion with NS boundary conditions on both cycles. We will 
do this in subsection 3.1 below, after reviewing the transformation from the flat torus to the 
singular surface £2,2- 

Alternatively, Refs. [9, 10] compute the reduced density matrix directly in terms of the 
two-point functions of Dirac fermions. It is clear that this calculation is for the ungauged 
fermion. The Hilbert space is factorized into left- and right-moving excitations, which is not 
possible for the Z2 gauged fermion dual to the boson; for example, the modular invariant 
partition function (which we will review below) does not factorize into contributions from 
left- and right-movers. In addition, the "modular Hamiltonian" , whose exponential forms 
the reduced density matrix (see equation (28) of [9]), consists of products of chiral fermion 
operators at different points. In the gauged theory, however, a single chiral fermion is not 
in the spectrum of local operators, as it is odd under Z2. Furthermore, the entanglement 
entropies calculated directly in this approach match those calculated via the replica trick. 

The free Dirac fermion is a consistent quantum theory, but it is not modular invariant, 
and so cannot be dual to the modular-invariant theory of the free boson at any radius. This 
explains the discrepancy between the Renyi entropies. However, there is a curious equality 
between the Renyi entropies for two intervals for the free Dirac fermion and the self-dual 
boson. We now turn to explaining this fact. 

3. Renyi entropies for arbitrary intervals 

The strange coincidence in the entanglement Renyi entropies for two intervals (N = 2) and 
all n, between the free Dirac fermion and the self-dual boson, challenges the ability of the 
entanglement entropies to distinguish theories. Before making any sweeping claims, we should 
compare the values of S n (A) for N > 2 in the two theories. Agreement or disagreement 
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of S n (A) between the two theories amounts to agreement or disagreement of the partition 
functions Z n ^ on the singular Riemann surface £ n jv- We will find on quite general grounds 
that the theories agree for n = 2 and any N (a new result), as well as for N = 2 and any n 
[8, 9, 10], for reasons which fail when n > 3 and N > 3. When n = N = 3 we check by direct 
numerical computation that the Renyi entropies differ. Hence the full set of Renyi entropies 
does distinguish between the theories. Along the way we will discover a surprising relationship 
between the theories, which will explain in a simple way why certain Renyi entropies agree. 

To compute the partition function, we can use a Weyl transformation to map the metric 
ds 2 on our singular surface T, n ^ to a non-singular fiducial metric ds 2 = e~^ds 2 [16, 15, 7, 6]. 
In this case, 

Z n ,N = e SL Z UtN . (3.1) 

Here Z is the partition function of the CFT with the fiducial metric, and Sl is the Liouville 
action 



c 



S l = ^JV9 [f'da^t + 2R^ . (3.2) 

The Liouville action depends on the CFT solely via its central charge. Since the Dirac fermion 
and free boson both have c = 1 , agreement of the Renyi entropies is equivalent to agreement of 
the partition functions on the non-singular Riemann surface. For the fermionic theory, it will 
also be important to keep track of boundary conditions around the various non-contractible 
cycles. 

In Appendix C, we apply the technology developed in this section to the compact boson 
at arbitrary radius, giving a relatively simple derivation of Calabrese, Cardy, and Tonni's 
result (2.4) for its Renyi entropies. 

3.1 Torus partition functions and momentum lattices 

The Riemann surface S^n has genus g = (N — l)(n — 1). We begin with the simplest non- 
trivial case, £2,2) which is a torus. The modular parameter r of the torus depends on the 
cross-ratio x of the endpoints of the intervals (the branch points in £2,2); defined in (2.5), 
which lies in the range < x < 1 when both cuts are on the real line. The relation between 
x and r is 

%F(l-x) . . /II \ 2 . , 

F(x) = 2 F 1 ( -,-,l;x) = -K(x), 

(3.3) 



F(x) ' y J z 1 \2' 2' ' J n 



where 2F1 is the usual hypergeometric function and K{x) is the complete hyperelliptic inte- 
gral. The expressions (3.3) are actually valid for arbitrary complex x, but for < x < 1, r 
is imaginary and the torus is rectangular. In this case, the partition function on a flat torus 
depends only on the spectrum of scaling dimensions A of the CFT: 

Z 2j2 = tr e - 27TT2H = e 2nT2 ^ 12 e~ 27TT2Am . (3.4) 
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The cycle playing the role of the spatial circle here is the one that, on £2,25 encircles one of 
the cuts while staying on one sheet, while the one playing the role of the time circle encircles 
the two middle branch points, passing from one sheet to the other. Based on the boundary 
conditions explained in subsection 2.5, both cycles have periodic boundary conditions for the 
fermions on £2,2, which corresponds to antiperiodic (NS) boundary conditions on the flat 
torus. Hence the trace is over NS-NS sector states and does not include a factor of (— 1) F . 

From the agreement of 82(A) for general x between the fermion and the self-dual boson, 
it follows that the two theories have the same partition function on any rectangular torus, 
and hence that they have identical spectra of scaling dimensions. In this subsection we will 
explain this agreement; in the rest of this section we will then use what we've learned to 
explain the agreement for other values of n (with N = 2), and find out to what extent it 
generalizes to other values of (n, N). 

Let us first recall the structure of the operators in the two theories. The Dirac fermion 
consists of separate left- and right-moving Weyl fermions 4>l,r, which have conformal weights 
h = 1/2 and h = 1/2 respectively. The general operator is a product of distinct operators 
of the form d n ipL, d n ^L, d n ipR, B n $R, where n = 0, 1, — The theory has conserved left- and 
right-moving fermion number currents. 

The compact boson can be split into left- and right-moving bosons Xl,r- The exponential 
operators, which create winding and momentum ground states, are of the form e lk LX L +ik R x R 3 
The left- and right-moving momenta Ul,r are not independent, but are elements of a joint 
lattice r( b ), which at the self-dual radius is as follows: 



The exponential operator has conformal weights h = &f,/2, h = fcJj/2. The total momentum 
and the winding number are given in terms of Ul,r by n,w = [k^ ± fcjj)/\/2. The general 
operator is a product of an exponential and derivative operators d n Xi, 8 n XR, n = 1,2,.... 
(The self-dual boson actually has a larger, SU(2) x SU(2) symmetry group, but we will have 
occasion only to use its momentum and winding [7(1 ) x U(l) subgroup.) 

For low-lying operators, it is straightforward to see by inspection that the spectra of 
the two theories are the same. For example, both theories have 4 dimension-1/2 opera- 
tors (fermion: ipL,R, $l,r', boson: the (n,w) = (±1, 0), (0, ±1) exponential operators) and 6 
dimension-1 operators (fermion: tpLtpR, "ipL^R, ipL,R^L,R] boson: the (n,w) = (±1, ±1) expo- 
nentials and OXl, OXr). Furthermore, it is possible to match not only the scaling dimensions 
but also the two theories' respective U(l) x £7(1) charges, i.e. to establish a one-to-one corre- 
spondence such that the left- and right-moving fermion numbers match n and w respectively. 
On the other hand, the spins definitely cannot be made to match, given that one theory 
contains fermions and the other doesn't. 

3 There are also cocycles, which we neglect here since they do not contribute to the scaling dimensions which 
are our main interest. In Section 4, where we will study correlators, we will include them. 




(3.5) 



- 13 - 



It is not necessarily clear from these low-lying examples, however, what the general 
pattern is. The mystery is readily solved by bosonizing the fermion. 4 In its bosonized form, 
the Dirac fermion consists of left- and right-moving bosons Hl^ r , related to the elementary 
fermionic fields by 

TpL = e iHL , Vi? = e iHR . (3.6) 

Just as for the self-dual boson, the general operator is written as a product of an exponential 
operator e lk LH L +k R H R anc j d er i va ti ve operators d n Hi,d n H R . But in this case the momenta 
^lr, which are the left- and right-moving fermion numbers, are independent integers; in 
other words, the momentum lattice is simply 

r (f) = Z 2 . (3.7) 

(Note that we use the superscript (f) to refer to the Dirac fermion theory even when we are 
working with its bosonized form.) The two lattices are shown in figure 4, and it is immediately 
seen that they are related by a 45° rotation. This rotation matches the U(l) x U(l) charges 
of the respective theories to each other: 5 

It also preserves the lengths of vectors defined with respect to the Euclidean inner product 

k ■ k' = k L k' L + k R k' R , (3.9) 

and therefore the scaling dimensions A = k ■ k/2 of exponential operators; however it does 
not preserve the Lorentzian inner product 

k o k' = k L k' L - k R k' R (3.10) 

which gives their spins s = kok/2. e In addition to having momentum lattices that are related 
by a rotation, we will show in Section 4 that the two theories have the same set of derivative 
operators. Together these two facts explain the matching of the spectra of scaling dimensions 
as well as U(l) x U(l) charges. 

Let us return to the torus partition function, which can be written in terms of a sum over 
the momentum lattice; this will be a useful warm-up for the higher-genus partition funcitons 



4 By "bosonizing" we mean representing the chiral fermion operators as exponentials of chiral boson oper- 
ators. The chiral bosons are denned via their OPEs. This is distinct from the bosonization of the modular- 
invariant fermion theory. 

5 We could just as well make other choices, like (k^p, ) = (n, —w) or (w, n), but these are all related by 
automorphisms of the two theories, and therefore equivalent. 

6 The Lorentzian inner product is perhaps more familiar in the context of momentum lattices. It is with 
respect to this inner product that, in string theory (for example in a Narain compactification), one requires the 
lattices to be integral (for mutual locality of operators), self-dual (for modular invariance under r — > — 1/t), 
and even (for modular invariance under r — > r + 1). and are both integral and self-dual, but only 
r (b) is even. 



- 14 - 



Dirac fermion Self-dual boson 




D I 2 -2-1012 



Figure 4: The momentum lattices for the ungauged Dirac fermion (left) and for the self-dual boson 
(right). We have bolded the canonical generating vectors for both lattices. Notice that they are related 
by a 45° rotation. 



we will consider in the rest of this section. We will not assume that r is imaginary. For the 
fermion we have 

Z^lir, f) = exp (iiTTik ok — ixr^k ■ k) . (3-H) 

The factor of | ^(-t - ) | 2 accounts for the sum over all the possible derivative operators that can 
multiply a given exponential (i.e. the oscillators, in terms of states). For the boson we have 
almost the same formula: 

( r , f) = |2 exp (iirTik ok — 7TT2 A; • k) . (3-12) 

^ T '\ fcer(t) 

Again, since r( f ) and r( b ) are related by a transformation that preserves the Euclidean inner 
product (3.9), the two partition functions will agree precisely when t\ = 0. In fact, this will 
work for any two lattices that are related by such an orthogonal transformation. However, if 
one restricts to integral self-dual lattices in two dimensions, then r( f ) and r( b ) are the only 
examples related in this way. 

From (3.11) we can easily recover the CFH result (2.3) for N = n = 2, as promised in 
subsection 2.5. Taking r imaginary, and using (3.3), we have 

Zvfr) - V ( T )2 - v^(i^y ) ■ (3 - 13) 

The Liouville action (3.2) for the Weyl transformation from the flat torus to the singular 
surface S2,2 was computed by Lunin and Mathur [16]: 

e^(M^)- 1/4 (^) 1/12 . 
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We obtain 



yielding (2. 3). 7 

In order to go to n > 2 and/or N > 2, we need to consider the partition functions of the 
theories on higher-genus Riemann surfaces, which we will do in subsection 3.3. To have the 
necessary language, however, we first need to review some algebraic geometry. 

3.2 Some algebraic geometry background 

In order to set up the computation of the partition functions Z n ^ n , we will describing the 
particular Riemann surfaces we are studying, and reviewing some basic facts about Riemann 
surfaces that we will need. More complete reviews of the relevant mathematics can be found 
in [11, 12, 17]. 

3rd sheet 




2nd sheet 




-<-n 



1st sheet 




A-cycles B-cycles 
■I.-.-.-H- -b* 



Figure 5: The Riemann surface £3,3 with a canonical basis of A- and B-cycles. The solid black lines 
are the branch cuts, oriented so that approaching from below takes one up a sheet and approaching 
from below takes one down a sheet. The blue dashed lines are A-cycles and the green solid lines are 
B-cycles, with notches corresponding to the index (eg 03, 04 lie on the second sheet). We have pulled 
the B-cycles off of the branch points for clarity. 



7 We can also calculate the result for the modular-invariant gauged Dirac theory. Using the fact that, in 
that case, Z = ' ?3+ 2 '^ +1?4 , we find J~2(x) = g ; which agrees with (2.4) for n = 2, 77 = 2, using the 

q2/V\ i „q2/^\ i „q2/ 



resummation identity #|(t) + 0|(t) + i?|(r) = 20 3 (2t)0 3 (t/2). 
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The singular Riemann surfaces were described in Section 2.4, and are illustrated in Fig. 
5 for n = N = 3. The Riemann surface can be represented by the algebraic curve [18, 19] 
(see also Appendix C of [6]): 

N-i 

y n = II ( z - u ^ z - v k) n ~\z - u N ) (3.16) 
k=l 

where Uk,Vk all lie on the real line. The kth. branch cut lies between branch points atUk,Vk £l 
with Uk < ffc. The residual conformal invariance can be used to fix the location of three of 
those points. For instance, one can set the first interval to lie between u% = 0, V\ = 1, and 
the final interval to end at vn = oo. 

The Riemann-Hurwitz formula gives the genus of this Riemann surface as g = (N— \){n — 
1). On any Riemann surface, one can write down a canonical basis of "A-cycles" ai=i... g and 
"B-cycles" &i=i ; ... ff with intersection pairing aj • bj = 5ij = —bj ■ en. For the curve (3.16), 
the canonical basis we will use is as follows: labeling the cycles by sheet s = 1, . . . ,n — 1 
and cut k = 1,...,N — 1, cycles cik+(s-i)(N-i) encircling the first k cycles on sheet s in a 
counterclockwise direction; cycles bfc+(s-i)(Ar-i) pass from to u^+i on sheet s, and return 
on sheet N [18, 19]. This basis is shown in figure 5 for the case n = N = 3. Based on 
the boundary conditions for the Dirac fermion explained in subsection 2.5, all of these basis 
cycles carry NS boundary conditions, which we can see as follows. First, each one crosses 
an even number of sign flips, and hence carries periodic boundary conditions on this singular 
surface. Second, each one has odd winding number. Hence, under a Weyl transformation that 
takes the surface to a non-singular one, the fermion will be antiperiodic along the geodesic 
representative of each cycle. 

There is a canonical basis of holomorphic differentials Q such that 



Cj = Sij ; f (j = SUj (3-17) 

a, J hi 

Qij is known as the period matrix; it is symmetric, and the imaginary part of Qij is a positive 
definite quadratic form. 

A natural but noncanonical basis of holomorphic differentials is [18, 19]: 

oj j+sm = - dz, j = l,...,N-l, 8 = l,...,n-l, (3.18) 

where y is given in (3.16). The period integrals are: 

Aij = & u)j, Bij = <L Uj, = A~ 1 B. (3.19) 

J ai J hi 

The canonical basis is then Q = ^^^k^kj ■ ^ s a warmu P) one can use the basis (3.18) to 
reproduce (3.3) for genus one, where fin = r. 
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3.3 Partition functions at higher genus 

The partition function of the self-dual boson can be computed following [11, 12, 20]: 



£00 



7T 



2Vol(£)Im det 



det' V 2 







(0|2«) 



(3.20) 



Here Vol(S) is the volume of the surface; e* G Z^ means e* € {0, 1/2}; det' V 2 is a suitably 
regularized determinant of the Laplacian; the factor of \/2ir comes from the path integral 
over the zero mode of the boson; and $ is the Riemann-Siegel theta function: 



r) 



r*2 



= ^ exp ^(n + e^ttijin + eiy +2iri(n + ei) k (z + e 2 ) k 



(3.21) 



The Dirac fermion with antiperiodic (NS) boundary conditions about all basis cycles 
can also be computed following [11, 12]. Those works do the calculation in the R = 1/2 
presentation of the dual boson, which can be written as a sum over spin structures of the 
fermion partition function. One can extract from this the partition function in the purely NS 
sector: 

r o 



29Vol(S)Im det Q 
det' V 2 



i) 







(o|n) 



(3.22) 



(which, again, includes the integral over the zero mode of the dual boson, equal to ir). 

To compare these partition functions, we first write = A + iK,, with A, K symmetric 
and real and K positive-definite. Using the results in Appendix A, we can rewrite (3.20) as 
a sum over g copies of the momentum lattice r( b ) for the self-dual boson, fcx, db /cr E \/2Z: 



^(b) 



\ l 2 V°K^)I™^ de t ^ eX p {jiiAijk 1 o k j - nK^k 1 • k j ] 



det' V 2 



(3.23) 



fce(r( b ))s 



Similarly, the expression (3.22) can be written directly as a sum over g copies of the momentum 
lattice r( f ) = Z 2 , 



2(f) 



7T 



2ffVol(S)Im det fi 



det' V 2 



y~] exp (niAijk 1 o k? - -KK^k 1 ■ k j ) . 
fce(r( f ))s 



(3.24) 



Since r( b ) is related to r( f ) by an orthogonal rotation, (3.23) and (3.24) will agree (up to a 
factor of 2( s_1 )/ 2 ) when Am = — that is, when the period matrix is imaginary. The factor 
of 2(f~ 1 )/ 2 can be absorbed in a constant shift of the dilaton coupling 



5S, 



dilaton 



(3.25) 



where R^ is the worldsheet curvature. Furthermore, with the result as stated, this factor 
will only contribute to the non-universal part of the Renyi entropy as a constant term, and 
will cancel out of the mutual Renyi information entirely. 



- 18 - 



The equality of the two partition functions for ReO = is the central result of this 
subsection. We will now return to the Riemann surfaces T, n ^ involved in the calculation of 
Renyi entropies, and ask whether £1 is in fact imaginary for them. 

3.4 Symmetry constraints on the Renyi entropies 

For the Riemann surfaces £ ni Ar , we will show that the real part of the period matrix vanishes 
identically when N = 2, and also when n = 2, but not otherwise. Thus all of the Renyi 
entropies for two intervals will agree between the free Dirac fermion and the self-dual boson, 
consistent with the results in [10]. In addition, the second Renyi for any number of intervals 
will also agree; in other words, the CFH result (2.3) for the Dirac fermion also applies to the 
self-dual boson for n = 2, which we believe is a new result. 

These facts arise because of the symmetries of the underlying Riemann surface. The 
Riemann surface S nj 7v is a cyclic branched cover with all of the branch points on the real 
line. Thus cyclic permutations of the sheets of the branched cover, as well as with complex 
conjugation, are symmetries of the surface. Together these symmetries form the dihedral 
group D n . The period matrix transforms in a specific reducible representation Rq of D n . 

These symmetries map A-cycles to ^4-cycles and .B-cycles to S-cycles. Thus, they preserve 
the canonical structure of the period matrix. The upshot is that the period matrix should 
satisfy the conditions 



where g is the (N — l)(n — l)-dimensional representation of the cyclic permutations acting on 
the 6-cycles, and r is the representation of the complex conjugation operators. This implies 
that: 



In other words, every nonzero component of A should lie in the one-dimensional irreducible 
representation of D n with g acting trivially and r = — 1; every nonzero component of K 
should lie in the trivial representation of D n . Our goal is to deduce the number of times these 
irreducible representations appear in Rq. 

To begin with we wish to deduce the representation of D n on the 6-cycles, decomposed 
into irreducible representations of D n , which we review in Appendix B. Consider first the 
case of 2 cuts. The generator g of the cyclic group Z n simply permutes the sheets of the 



8 General treatments of real algebraic curves can be found in [21, 22]. These use a different canonical 
homology basis than the one discussed here, related by a symplectic modular transformation. Since the free 
Dirac fermion theory is not invariant under modular transformations, the partition function in that basis will 
not directly give the related Renyi entropy. 



(3.26) 



g T Kg = K, tKt = K, g T Ag = A, tAt = -A. 



(3.27) 
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branch cover. This, it acts on the -B-cycles b, 



s=l 



/-I -1 

1 
1 



\ 



i-i as 

-1 -1\ 





(3.28) 



1 J 



The complex conjugation t takes the complex conjugate of the first sheet; for n even, it 
also acts as complex conjugation on the + l)st sheet. Acting on the additional sheets, it 
exchanges the kth sheet with the complex conjugate of the (n — k + 2)nd sheet. We can thus 
write the action on the 6-cycles as: 
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The eigenvectors of g with eigenvalue A have the form 

A n-3 



A 

V i / 



such that 



k=l 







(3.29) 



(3.30) 



(3.31) 



For n odd, the solutions are A& = e 2mk / n . Based on the representations listed in Appendix 
B, we deduce: 

R n = R 1 ®R 2 ®...®R n=1 (nodd). (3.32) 

2 

For n even, the solutions are A& = e 2nlk / n . These include pairs of roots of unity and 
= — 1- F° r ^ = — 1) t ne corresponding eigenvector of g, v\_ 1} is an eigenvector of r with 
eigenvalue 1. Thus, in this case 

Rn = S © Ri © R 2 . 



R n~2 
2 



(n even). 



(3.33) 



In the case of tr p n for two cuts, the period matrix transforms as the symmetrized product 
(Rn®Rn) S ym- To work this out, we need to decompose products of irreducible representations 
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of D n . This can be simply worked out from the specific matrix forms of g, r given in Appendix 
B: 



Rj ® Rk = Rj+k © Rj-k m 



Rj ® Rj = R 2j © 1 © T 



Rj <g> S = R n-2k 



2 



S®S = 1 



(3.34) 



The first two hold for n even or odd. The representation T only occurs in the second line. 
However, it corresponds to the antisymmetrized tensor product of Rj with itself, and does 
not appear in (Rq <g) Rn)sym- Therefore, the real component of the period matrix is zero in 
the two-cut case. The Renyi entropies of the free Dirac fermion for two disjoint intervals are 
then identical to those for the self-dual boson, consistent with the results of [10]. 

Next, let us consider the case of 3 or more cuts. As discussed in subsection 3.2 and 
Fig. 5, the basis of -B-cycles can be indexed as 6^. + ( s _ 1 wjy_ 1 ). In this case D n acts as the 
representation Rq on the indices s, and the full representation is (B^Z-^Rn- The period 
matrix can be written in (n — 1) x (n — 1) blocks. The N — 1 diagonal blocks are each in 
the representation (Rq ® Rn)sym, and so have no real parts. The symmetry of £1 equates the 
(A — 1)(A — 2)/2 off-diagonal blocks below the diagonal to those above the diagonal. Each of 
these blocks transforms in the unsymmetrized representation Rq £§> ifo. There are (n — 2)/2 
occurrences of the irrep T in each off-diagonal block when n is even, and (n — 1)/2 occurrences 
in each off-diagonal block when n is odd. Thus, the expected number of real components of 
the period matrix is: 



The second line implies a new result: for the second Renyi entropy 82(A), the CFH result for 
the fermion (2.3) [8, 9, 10] applies also to the self-dual boson. 9 

Outside of the two families A = 2, any n and n = 2, any N, symmetry constraints 
do not prevent a real part of the period matrix. We have computed the period matrix and 
partition functions for the case A = 3, n = 3 using Mathematica, both by using the built- 
in function SiegelTheta and by explicitly computing the theta functions, truncating the 
infinite sum (3.21) at rij = 5. The two methods agreed with each other to six decimal places. 
For a variety of locations of the branch cuts of the Riemann surface, we found a single real 
component, as predicted in (3.35), which varied with the locations of the branch points, and 
we found that the partition functions of the free Dirac fermion and self-dual boson differed 

9 A very similar calculation gives us the number of imaginary components of the period matrix. These must 
lie in the trivial representation 1 of D n . For n odd, this gives N[N — l)(n — l)/2 imaginary components for 
71 odd, and N(N — l)n/4 components for n even. 



(A - 1)(A- 2)(n- 1) 



components for n odd 



4 



(N — l)(JV-2)(n-2) 



components for n even 



(3.35) 
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at the percent level. As a check, if we set the real component of the period matrix to zero 
and inserted this into the partition functions (3.20,3.22), we found that they agreed to six 
decimal places. We conclude that, in the end, the entanglement spectrum does discriminate 
between different theories. A natural question to ask now is, what structures in the CFT are 
being singled out by the n, N > 3 Renyi entropies? We will give a conjecture at the end of 
the next section. 



4. A real duality 

We saw in the last section that the free Dirac fermion theory and the self-dual compact boson 
have identical spectra of scaling dimensions. We also saw that this coincidence can most 
easily be understood by bosonizing the fermion to obtain left- and right-moving bosons with 
a momentum lattice, Z 2 , that is related to that of the self-dual boson by a 45° rotation 

m mx f k {h) + k {h) k {h) -k {h) \ 

k?,4 } ) = ( , (4-1) 

which preserves the scaling dimension A = {k\ + kj^/2, but not the spin s = (k\ — kj^/2. 

Unlike a true duality (such as T-duality), which acts separately on the left- and right- 
movers of a CFT, the 45° rotation we are discussing mixes left- and right-movers. Nonetheless, 
it does give identical results for the partition functions on whole families of Riemann surfaces. 
Friedan and Shenker [23] have suggested that the partition function, understood as the section 
of a line bundle over the moduli space of Riemann surfaces, for all genera, should define the 
conformal field theory (via its factorization limits). In this spirit, it is worth asking whether 
the 45° rotation (4.1) might preserve some of the structure of the two CFTs, such as correlators 
and OPEs. The purpose of this section is to show that a large amount of the structure — in 
fact, all correlators and OPEs on the real line — is preserved. 

Consider for example the two-point function of elementary fermionic fields ipLi^L = 
e ±iH L . 

Z\ Z2 

These operators are mapped by (4.1) to the operators of the self-dual boson with momenta 
±1 and no winding; these are spinless operators, whose two-point function is 

e i(X i ( 2l )+X fl (2 1 ))/ v / 2 e -i(X i ( Z2 )+X ii (5 2 ))/v / 2^ _ 1 _ (43) 



\Zl - z 2 



While the correlators do not match for general positions, we see that they do match whenever 
z \ — Z2 is real and positive. Of course, two-point functions of primary operators are dictated by 
their scaling dimensions, so the nontrivial question is whether the agreement extends to non- 
primaries and to higher-point functions. In fact, as we will show, it does. The only restriction 
is that they must be placed on the real axis (or any common horizontal line) in the same order 
that they are written in the expectation value; that is, for the correlator (Oi(zi)02(z2) ■ ■ ■) 
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we require the Zi to be real and z\ > Zi > • • • . We will prove this correspondence, which we 
call a "real duality", in subsection 4.1. We will begin with the exponential operators, which 
are related by (4.1), then explain how to generalize the correspondence to other operators in 
such a way that OPEs on the real line are preserved; this immediately implies that all n-point 
functions on the real line are preserved. 

In subsection 4.2, we will show that the component T xx = T + T of the stress tensor is 
preserved by this correspondence. Using the OPE, this implies that action of the mixed Vira- 
soro generators L n + L n are preserved, generalizing the matching of scaling dimensions. From 
this we learn what kinds of conformal transformations commute with the correspondence. 

Finally, in subsection 4.3 we will use the Schottky construction, which describes an arbi- 
trary Riemann surface as a quotient of the Riemann sphere by a discrete subgroup of SL(2, C), 
to explain how the real duality is related to the equality of partition functions between the 
two theories on Riemann surfaces with imaginary period matrices. 

Throughout this section we will refer to the ungauged Dirac fermion theory as the 
"fermion theory" (although we will mainly work in its bosonized form), and the theory of the 
self-dual boson as the "boson theory". We think our meaning will be clear. 



4.1 Correlators 

In this subsection we will define a one-to-one correspondence between the operators of the 
fermion and the boson theories, and show that under this correspondence arbitrary correlators 
in the two theories agree, as long as the operators are placed on the real axis (in the same 
order that they are written in the expectation value). To prove this, we will show that the 
OPEs on the real axis agree; the statement about correlators follows since an arbitrary re- 
point function can be reduced to 1-point functions by repeated application of the OPE. Since 
the OPE proof is a bit formal, it is perhaps useful to see the real duality in action first. So we 
begin in 4.1.1 by proving it by explicit calculation for exponential operators. Then in 4.1.2 
we explain how to generalize the one-to-one correspondence from exponential operators to 
general operators. Finally, we give the OPE proof in 4.1.3. 



4.1.1 Exponential operators 

In both the fermion and the boson theories, the exponential vertex operators include a co- 
cycle, which is necessary because the pure exponential operators have the wrong statistics. 
For example, in the fermionic theory e lHL and e lHR commute even though they represent 
fermionic operators, and in the bosonic theory e l ^ XL ^ +Xa ^ 1 ^/ y ^ and e % ^ XL ^ Z2 ^~ XR ^ Z2 ^ / ^ 
anticommute even though they represent bosonic operators. In the fermion theory, the com- 



- 23 - 



plete vertex operators are 10 

y( f ) — (_i^k R p Le ik L H L +ik R H R ^ (4-4) 

where pl is the left-moving momentum operator (the operator whose eigenvalue is ki), while 
in the boson theory they are 

y( b ) _ ^_-^{k L -k R ){p L +p R ) /2 e ik L X L +ik R X R ^ ^ 

We will now compute the correlators of these vertex operators, and prove the agreement 
claimed above. In the fermion theory we have 



x ' [ , otherwise 

(4-6) 

The sign factors (— l) k L k R come from moving the cocyles past the exponentials to the left 
until they hit the vacuum. The rest comes from the expectation value of the exponentials. 
In the boson theory the result is similar: 



(v i k \\z 1 ,z 1 )V i 



(b) 

fc 2 1 



Z2,Z 2 



, otherwise 



The momentum conservation conditions are linear, and therefore preserved by the rotation 
(4.1). The sign from the cocycles is clearly preserved by (4.1). Finally, when the Zj — Zj are 
all real and positive, the rest of the multiplicands collapse to [z% — Zj) k *' kJ , and the exponents 
are again equal under the rotation. So indeed (4.6) and (4.7) are equal under (4.1). 

4.1.2 General operators 

We now wish to extend the correspondence from exponential operators to general operators, 
which are products of the operators Vfc and derivative operators d n Hi, d n Hn (for the fermion) 
or d n XL , d n Xn (for the boson) . Operators that include derivatives are degenerate for fixed 
values of the momenta and scaling dimension, so we cannot be guided by matching those 
quantum numbers alone. Our guiding principle for dealing with such operators will be the 
following. Since the correlators of exponential operators match only when they are placed on 
the real axis (or any other common horizontal line, but for simplicity we will take the real 



10 With respect to the circle product k ok' = kLk' L — k R k' R , the lattice Z 2 for the fermion theory is integral 
and self-dual but not even (i.e. the theory contains fermionic operators). The standard cocyle prescription, 
described in textbooks, applies only to even lattices, and as far as we know there is no simple general expression 
for the cocycles for non-even lattices (see [24] for a discussion of this point). That this particular theory admits 
a simple expression for the cocycle is presumably due to the fact that its lattice, while not even, is simply 
related to an even one. 
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axis), we will work entirely on the real axis. Now, by the equations of motion for the fields, 
the holomorphic and antiholomorphic derivatives d, d can be replaced by the x-derivative 
d x = d + B: 

8 n H L = 3 n x H L , B n H R = d n x li R , (4.8) 

and similarly for Xl,X r . Having written all derivatives in terms of d x , we simply apply the 
same rotation in field space that yielded the map between the exponential operators, namely 

H L v>±={X L + X R ), H R <*±={X L -X R ). (4.9) 

Thus for example we have 

d n H L = d n x H L o ^d n x (X L + X R ) = ^(0 n X L + B n X R ) , 

d n H R = d n x H R O \%{X L - X R ) = ^{d n X L - d n X R ) . (4.10) 

It is useful to note that the rotation (4.1) on the momentum lattice, together with the 
rotation (4.10) on the derivative operators, clearly preserves the Zamolodchikov metric on 
the space of operators. Since, by definition, the Zamolodchikov metric is the correlator 
Gmn = (^m.( oo )-^n(0)) (where A' m is the operator A m in the z' = l/z frame), this implies 
that n-point functions where one of the operators is at infinity are also preserved by the real 
duality. We will also make use of the matching of the Zamolodchikov metric in subsection 
4.3. 



4.1.3 OPEs 

We will now show that the OPE of arbitrary operators at real positions z\ , Zi , with z\ > Z2 , 
is preserved, i.e. if <-> J< b ),£W <-> then 

^ i \z 1 ,z 1 )g^(z 2 ,z 2 )^^(z 1 ,z 1 )g^(z 2 ,z 2 ). (4.11) 

As mentioned above, the agreement of arbitrary correlators on the real axis follows from the 
agreement of OPEs. 

All of these composite operators are defined via normal-ordering; it will be useful to 
indicate this explicitly. Recall that in a free theory the OPE of operators : J-(zi,z~i) :, 
: Q(z 2 , z 2 ) : is derived by adding to : T{z\, zi)Q(z 2 , z 2 ) : all possible cross-contractions between 
J- and G, and then Taylor-expanding with respect to Z\ — z 2 and z\ — z 2 inside the normal- 
ordered product. We will first show that the cross-contractions match, then that the Taylor 
expansions match, when z\ — z 2 is real and positive 

In the fermion theory, a cross-contraction consists of replacing an Hl in J- and an Hi in 
Q with the propagator — ln(zi — Z2), or a pair of H R s with — ln(zi — z 2 ); similarly in the boson 
theory with Xl and X R . Let us first see how this works for derivative operators, then we 
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will give the general proof. For example, the OPE of the left-moving fermion number current 
8Hl with itself is 



: dH L (z x )dH L {z 2 ) : - - — ^ . (4.12) 



: dH L {z x ) :: 8H L (z 2 ) : = : dH L {z x )dH L {z 2 ) : - d 1 d 2 hi{z l - z 2 ) 

1 

{z\ - z 2 y 

The corresponding operator in the boson theory is the momentum current (dXi + BX R ) / y/2. 
Multiplying it by itself, we have four terms; two of them have mixed holomorphic and anti- 
holomorphic parts and therefore no cross-contractions, while the other two are of the same 
form as (4.12): 

i: 8X L {zi) + BX R { Zl ) :: dX L {z 2 ) + BX R (z 2 ) : 

= \: (dX L ( Zl ) + BX R { Zl )) (8X L (z 2 ) + BX R (z 2 )) : - 1 - 1 _ . (4.13) 
z z{zi - z 2 y z(zi — Z2J 

Comparing the right-hand sides of (4.12) and (4.13), clearly the operator parts match under 
(4.10), while their c- number parts are equal whenever the Zi are real. More generally, a 
contraction of d ni HL(z\) with d n " z K^z 2 ) gives 

- W H Zl - Z2 ) = (-^r^ + ^ -iy. 

while the contraction of the corresponding operators in the boson theory gives 

Vfl-lnf* z) 1 B^\n(z f (-l) ni (ni + n 2 -l)! (m + n 2 - 1)! 

--^ 9 2 ln(zx - z 2 ) - d 2 ln( Zl - z 2 ) - 2(zi _ Z2)ni+n2 + 2(z - 1 _ - Z2) n 1+ n 2 ■ 

(4.15) 

Again, these are equal when the zi are real. Clearly the same thing will hold for the contraction 
of B ni H R ( Zl ) with B n2 H R (z 2 ). 

For the proof that the OPEs of general operators match, we now apply the general formula 
for the cross-contractions in a free field theory (see for example equation (2.2.10) in [25]). In 
the fermion theory this is 

: ?{z\,z\) :: £(2:2,^2) : 



exp 



I^ 222 { HZl ~ ^M^SH^zV) + ^ - ^5H 5 R ( Zl )6H 6 R (-z 2 ))_ 

x :F( Zl ,z 1 )Q(z 2 ,z 2 ):, (4.16) 



Since the operators involved are contained entirely on the real axis (including all derivatives, 
when written using d x ), we can replace the integrals and functional derivatives with respect 
to z, z with ones with respect to x: 



: F(xi) :: Q(x 2 ) : ■ 



exp 



j d x ld x 2 H Xl -x 2 ) 



+ 



6H L ( Xl ) 5H L (x 2 ) 6H R ( Xl ) 5H r (x 2 ) 



: F(xi)Q{x 2 ) : . 

(4.17) 
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Similarly, in the boson theory we have 



: F{x{) :: Q{x 2 ) : : 



cxp 



dx\dx 2 ln(xi — X2) ( — r — r + 



5X L ( Xl ) 5X L (x 2 ) 5X R ( Xl ) 6X R (x 2 ) 



F(xi)G(x 2 ) : . 

(4.18) 



These clearly map to each other under (4.9). Note that we also need x\ > x 2 , otherwise the 
branch cut in the logarithm can lead to a mismatch (as between (4.2) and (4.3)). Finally, 
we note that moving the cocycle for Q through T gives a factor of (— l) k L k R in the fermion 
theory and of (—i)™" 7 ^ 6 [ n the boson theory, but these are equal under the correspondence. 

Finally, continuing to work on the real axis and to express derivatives using d x , it is clear 
that Taylor-expanding the normal-ordered operators obtained from the cross-contractions 
with respect to x\ — x 2 will commute with the map (4.9). Hence the correspondence we have 
described preserves the full OPEs of arbitrary operators on the real axis. 

4.2 Stress tensor and real conformal transformations 

The components of the stress tensor in the two theories are 

T (f) = A d H L dH L , f « = ~BH r BH r (4.19) 

and 

T (b) = _ l -dX L dX L , f ^ = - l -dX R dX R , (4.20) 

respectively. While these components do not individually map to each under the correspon- 
dence (4.10), their sum (which is the component T xx ) does: 

y(f) + f-(f) T {h) + T (b) . (4.21) 

The OPEs of T, T with other operators determine the action of the Virasoro generators: 

T(z)A(0) = z- (n+2) L n ■ A(0) , f(z)A(0) = £ z~^L n • A(0) , (4.22) 

n n 

where L n ■ A denotes the result of L n acting on the operator A via the state-operator map- 
ping (i.e. if the state-operator mapping maps \A) to A, then it maps L n \A) to L n ■ A). 
Since the OPEs between operators on the real axis are preserved by the correspondence, for 
corresponding operators A^ , ^4^ b ^ , we have 

J>-(«+ 2 ) (4 f ) + -^ f )(o) = (r«(x) + r( f )(x)) A^(o) 

n 

<-> (r ih \x) + f {h \x)^A ih) {o) 

= £ x~^ (4 b ) + L^) • ^ b ) (0) . (4.23) 
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Hence the action of L n + L n commutes with the correspondence. This statement generalizes 
the fact that the scaling dimension, which is the eigenvalue of Lq + Lq, is preserved. 

The actions of L n ,L n in turn determine how an arbitrary operator transforms under 
conformal transformations. Specifically, under an infinitesimal conformal transformation z \— > 
z + ev(z), we have: 

8A(z, z) = -eJ2 -j (v (n) (z)L n -i + (z)*L n -i) ■ A(z, z) , (4.24) 

n=0 

where v ^ = d n v. Since the action of L n + L n is preserved by the correspondence, if 
is real for all n (at the location of the operator) then 

<L4 (f) <L4 (b) . (4.25) 

Exponentiating an infinitesimal transformation such that v^ n \x) is real for all real x yields a 
finite transformation described by a real analytic function x'(x) with positive first derivative. 
Such "real conformal transformations" are compatible with the correspondence between the 
two theories, in the sense that if A^\x) o A^(x) then A'^(x') A^°\x'). Note that real 
conformal transformations also preserve the order of the positions of operators along the real 
axis. 

Just as the notion of a usual (complex) conformal transformation is local and can be used 
to do a coordinate transformation on a patch of a manifold, the same holds for real conformal 
transformations. A simple example is afforded by the cylinder defined by identifying the plane 
in the imaginary direction w ~ w + 27ri. The conformal transformation z = e w maps the 
cylinder to the plane with the origin removed. The map induces a real conformal transfor- 
mation from the real axis in the w-plane to the positive real axis in the z-plane. Therefore, 
since the correspondence between the fermion and boson theories preserves correlators on the 
real axis of the z-plane, it also preserves correlators on the real axis of the ^-cylinder. Such 
correlators represent, in the Lorentzian theory, equal-time correlators at finite temperature. 

4.3 Partition functions and the real duality 

In Section 3 we showed that the partition functions of the Dirac fermion and the self-dual 
boson agree on Riemann surfaces with imaginary period matrices. In this section we have 
shown so far that, under a certain one-to-one mapping between the operators of the theories, 
the correlators on the real axis agree. Both results are essentially due to the fact that the 
momentum lattices for the two theories are related by a rotation that preserves the dot 
product, k ■ k' = ki}i' L + kftk' R . In this section we will argue that these results can also be 
directly related to each other, since a genus-g partition function can be written in terms of 
2g-point functions on the plane. As we will review, there is a relationship between the period 
matrix and the positions of the operators which is such that the period matrix is imaginary 
if the positions are real. We conjecture that the reality of the potions is both sufficient and 
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necessary for the period matrix to be imaginary. If this conjecture is true, the agreement 
between the partition functions follows from the real duality. 

We first briefly review the relation between the genus-<7 partition function of a general 
CFT and 2g-point functions on the plane, closely following the discussion in the appendices of 
[26] . We begin with the Schottky construction, which describes an arbitrary genus-g Riemann 
surface as a quotient of the Riemann sphere by a discrete subgroup V of SL(2, C). This is a 
free group with generators ji, i = 1, . . . , g, that act as follows: 



7i(z) 



Pi 



(4.26) 



we have parametrized ji in terms of its attractive and repulsive fixed points Oj , r i and the 
dilatation parameter pi, which satisfies < \pi\ < 1. To obtain a fundamental domain for T, 
we remove from the plane, for each i, the discs 



Di 



where R; _,■ are chosen so that 



< R. 



R;R~ 



[Pi 



(4.27) 



(4.28) 



and so that none of the discs overlap. Consistency of these requirements places some restric- 
tions on the Schottky parameters a,i,ri,pi. To reconstruct the Riemann surface, one identifies 
the boundaries of by the action of 7$. There is a fundamental basis of A and -B-cycles, 
in which the A-cycles are represented by the dDi, and the B-cycles by lines connecting Di 
to D—i. With the restrictions mentioned above, the Schottky parameters cover the moduli 
space of genus-g Riemann surfaces (with some redundancies; for example, one can conjugate 
the entire Schottky group by an element of SL(2, C) without changing the Riemann surface). 

By the standard sewing construction, the partition function of a CFT on the surface 
obtained from the Schottky construction can be written as a sum of 2g-point functions on the 
plane. We will simply quote the result here; the detailed derivation can be found in appendix 
C of [26]. Let {^l m } be a basis of operators with conformal weights h m , h m respectively, and 
let G mn be the Zamolodchikov metric and G mn its inverse. Then 



mi,ni,...,m J ,fij \ i 



I ((ri - ai) L »{r t - ai) lo e L ^A mi {ri)) Un - a^in - a;) ifl 'e~ L ^ ^0*0 



(4.29) 



We will now show that (4.29) gives the same result for the partition function for the 
fermion and boson theories, whenever the pi,a{,ri are all real and the pi are all positive. 
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First, in that case the formula simplifies as follows: 



Z g {pi; at; n)= Hpf^G^isgn^ - a,)) 2 ^ 

mi,ni,...,m s ,n 9 \ i / 

x l^^n-a^+^e^+^AnM)) (in-a^+^e-^-^An^y (4.30) 

The factor of (sgn(rj — aj)) 2Sm ! arises from writing (ri — ai) L °(fi — ai) L ° = |r»— ai\ Lo+L ° sgn(rj — 
a ^L —L . s j nce anc [ j^ n ^ necessarily have the same spin (otherwise their Zamolodchikov 
inner product would be zero), and e ±( L i+ L i) doesn't change their statistics, we can write 
the sign factor as (sgn(rj — ai)) 2Sm i and pull it out of the expectation value. From the real 
duality — specifically, the agreement of scaling dimensions, the Zamolodchikov metric, the 
action of L n + L n , and correlators on the real line — we see that (4.30) gives the same result 
when applied to the fermion and boson theories. The factor of (sgn(rj — ai)) 2Sm i takes care 
of the fact that the A mi (ri) and A ni (a>i) are in the "wrong" order (for the real duality) when 
en > T{\ finally, the order of the multiplicands for different values of i inside the expectation 
value doesn't matter, since A mi and A ni are either both bosonic or both fermionic, so the full 
multiplicand is always bosonic. Since they have the same central charge and therefore the 
same Weyl anomaly, they will also have the same partition function on any surface related to 
this one by a Weyl transformation, including the constant-curvature one. 

We have now proven that the partition functions of the fermion and the boson theories 
are the same under two separate sets of conditions on the moduli of the Riemann surface: 

1. when the period matrix Q is imaginary (in Section 3); 

2. when the Schottky parameters ai,ri,pi are all real and the pi are positive (just above). 

While we are not aware of a proof in the mathematical literature, it seems very likely that 
these two sets of conditions are actually equivalent, i.e. that the period matrix is imaginary 
precisely under the conditions (2) on the Schottky parameters (or rather, when the pi are 
positive and the a%,ri can be chosen to be real using the SX(2, C) freedom). 11 The relation 
between the Schottky parameters and the period matrix (in the basis of A- and S-cycles 
described above) is known explicitly, but is somewhat complicated: 



(Oi- l(ai)){n -i(n)) 
(oi - l(ri)){ri - 7(a»)) 



e" " = Pi 



TT J r (Mi). (4-31) 

Here (7^) \ Y/(~jj) is the set of all 7 G T, written as words made from the letters Jk=i,...g an d 
their inverses, such that the first letter is not 7^ or 7 i ~ 1 and the last letter is not jj or lj l ■ 



L We would like to thank M. Gaberdiel, R. Volpato, and X. Yin for helpful discussion on this point. 
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From (4.31) we can almost prove one direction of the equivalence. If the Schottky parameters 
are real, then obviously e 27r ij is real for all If in addition pi > 0, then e 2mn " > 0, hence 
fiji is imaginary (the product over 7 is a perfect square, hence positive, since the multiplicand 
has the same value for 7 and 7" 1 ). Presumably the same can be shown for the off-diagonal 
components of £1. Showing the converse, that if £1 is imaginary then the pt are positive and 
the ai,ri can be chosen to be real, seems more challenging, and we will not attempt it here. 

5. Discussion 

5.1 Generalizations 

5.1.1 Other configurations and states 

So far in this paper we have taken the field theories being studied to be in their ground states. 
However, it is straightforward to generalize the analysis to finite-temperature states. In this 
case the Euclidean spacetime that gets replicated in the replica trick is periodically identified 
in the Euclidean time direction (with NS boundary conditions for fermions), giving a cylinder. 
Including the points at spatial infinity, this is a sphere, so the replicated surface has the same 
topology as at zero temperature. Although it has a different complex structure from the zero- 
temperature case, the same basis of cycles can be used, and the dihedral representation theory 
argument given in subsection 3.4 goes through as before. Hence, just as at zero temperature, 
this Riemann surface has an imaginary period matrix for N = 2 and for n = 2, and therefore 
that the Renyi entropies for the Dirac fermion and the self-dual boson agree in these cases. 

Another generalization is to quantize the theories on a circle, rather than a line. In 
this case, the replicated surface is periodic in the spatial direction. At zero temperature 
we again have a cylinder, and if we put NS boundary conditions on the fermion then the 
analysis of the previous paragraph shows that we will again get agreement for N = 2 and 
for n = 2. On the other hand, if we consider the theories on a circle at finite temperature, 
then the Euclidean spacetime is a torus. The ra-sheeted replicated surface now has an extra 
2n cycles, which transform in the fundamental representation of the dihedral group. The 
representation theory is therefore the same as if we had added another cut (plus one extra 
trivial representation). Hence the Renyis will agree for N = 1 and for n = 2. 

Returning to the theory on the line at zero temperature, a different generalization is to 
intervals in spacetime that do not lie on a constant-time line. These are more difficult to 
compute, since the usual replica trick cannot be applied. However, it is possible that these 
quantities are related, perhaps by some sort of analytic continuation, to Euclidean partition 
functions where the endpoints Uj,Uj of the intervals are moved off the real axis. Since our 
explanation of the agreement between the Renyis for the Dirac fermion and the self-dual 
boson crucially required those branch points to be real (in particular in our analysis of the 
symmetries of the relevant Riemann surface in subsection 3.4), it seems very likely that the 
non-equal-time Renyis will indeed distinguish between the theories, even for two intervals. 
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5.1.2 Other pairs of theories 

The free Dirac fermion and the compact boson are among the simplest quantum field theories 
one can study. An obvious question is whether the coincidences we have found — concerning 
entanglement entropies, partition functions, and correlation functions — can occur for more 
complicated theories, or whether they are in some sense artifacts of these theories' simplicity. 

One generalization to a class of more complicated theories follows straightforwardly from 
our analysis. The key relationship between the Dirac fermion and self-dual boson theories, 
which allowed us to show that their partition functions agreed for imaginary period matrices 
(hence their Renyi entropies for N = 2 and for n = 2), as well as to prove the real duality, was 
the fact that their momentum lattices T™\ are related by a transformation that preserves 
the Euclidean inner product k ■ k' = kik' L + knk' R . It is clear from our analysis that any two 
theories that can be described in terms of left- and right-moving bosons on momentum lattices 
related by an orthogonal transformation will enjoy the same set of coincidences. Trivial 
examples include T-duality (which takes kn — > — /cr) and parity (which exchanges k^ and 
kji), but more interesting examples will mix left- and the right-movers. For integral self-dual 
two-dimensional lattices, and r( b ) furnish the only such example, but presumably with 
more bosons there are more examples. In fact, it would be interesting to see whether there 
are pairs of even self-dual lattices related in this way. One could also consider lattices that 
are not integral or not self-dual. 

More generally, it would be interesting to study whether similar "partial" dualities can 
occur for theories that are not described in terms of free bosons, or in higher-dimensional 
theories. 

The special role of dynamics on a codimension-one surface has the flavor of boundary 
conformal field theory [27], of which a standard example is taking a CFT on the upper half 
plane with some boundary conditions on the real line. An exactly solvable BCFT in terms 
of an interacting self-dual scalar was described in [28], and was nontrivially fermionized in 
[29] (leading to a different theory than the non- modular invariant Dirac fermion). It would 
be very interesting to know if this nontrivial exact equivalence of theories is related in some 
way to our real duality, which is of course not an exact equivalence. 

5.2 Connections to larger issues 

Taking account of what we have learned in this paper, we return in this final subsection 
to the questions we posed at the beginning of the paper: Are entanglement entropies in 
quantum field theories (or, more precisely, their finite parts) universal quantities, and do they 
distinguish between theories? 

5.2.1 Position- vs. momentum-space entanglement 

The results of section 2 are consistent with the universality of the entanglement spectrum of 
the reduced density matrix for spatial subsets of a field theory. More precisely, the cutoff- 
independent quantities such as mutual informations are independent of the specific Lagrangian 
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presentation. This is consistent with the statement that a conformal field theory is defined 
by the spectrum of local operators and the operator product expansion. 

We have shown the equivalence of the real-space entanglement spectra between boson and 
fermionic presentations at a specific point in the moduli space of c = 1 conformal theories, 
where both theories are free. There is a line of conformal field theories, which corresponds to 
different radii of the free boson, and to a varying four-fermion coupling in the dual, modular- 
invariant fermionic theory (here we do not mean the ungauged Dirac fermion related by "real 
duality" ) . We expect that the entanglement spectra of this whole line of theories are invariant 
under the bosonization map. 

In contrast, one may consider the entanglement of regions in momentum space, as dis- 
cussed in [30]. This entanglement is an important aspect of Wilsonian renormalization, in 
which ultraviolet degrees of freedom are integrated out, or traced over; in any interacting 
theory, the ultraviolet and infrared degrees of freedom are entangled in the ground state, and 
the state of the IR theory is described by a density matrix. However, this entanglement is not 
universal in the same sense; rather it depends very much on one's choice of presentation. In 
the case of Bose- Fermi duality for arbitrary boson radius, if we choose the ultraviolet degrees 
of freedom to correspond to bosonic oscillators at high momentum, then the ground state 
will factorize between ultraviolet and infrared as the theory is free. However, if we choose 
the UV degrees of freedom to correspond to fermionic oscillators at high momentum, then 
the four-fermion interactions guarantee that the ground state will be highly entangled be- 
tween momentum scales, as is apparent by studying the explicit construction of the bosonic 
ground state in the interaction picture of the fermion theory [31]. In general, this presentation 
dependence is related to the fact that in calculations of real-space entanglement, the cutoff- 
dependent terms are scheme-dependent — they depend on the details of how one partitions 
the theory between IR and UV degrees of freedom. 

This is not to say that momentum-space entanglement is not useful; it is an important fact 
about integrating out UV degrees of freedom in an interacting theory [30], and is a measure of 
the interactions of a given set of degrees of freedom. But the real-space entanglement appear 
to be the right tool for characterizing theories in an invariant manner. 

5.2.2 Do Renyi entropies distinguish theories? 

Our study of the free fermion and self-dual boson has given evidence that entanglement 
entropies do indeed distinguish theories, if one includes regions with enough components. As 
we mentioned in Section 4, this question is related to an old program of Friedan and Shenker 
[23], who proposed that the set of partition functions on Riemann surfaces of all genera, as 
a function of the moduli of those Riemann surfaces, might completely characterize modular- 
invariant conformal field theories. The essential point is that the factorization limits build 
up these Riemann surfaces in terms of correlation functions on the sphere. The replica trick 
relates this proposal to the attempt to characterize conformal field theories via their Renyi 
entropies. The Renyi entropies, however, only depend on the Riemann surfaces on a slice 
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through the full moduli space. The question remains as to how much information about the 
conformal field theory can be extracted from this restricted class of partition functions. 

Finally, we would like to point out that there exists a large class of theories for which the 
distinguishing ability of entanglement entropies, along with the Friedan-Shenker program, 
fails badly in a specific limit. These are theories at large c and strong coupling — very far 
from the free c = 1 theories we've been considering so far. Specifically, any holographic CFT 
whose dual is Einstein gravity (possibly coupled to some matter) will have the same partition 
function at leading order in 1/c on a given Riemann surface. The reason is that the partition 
function is determined by the solution to the Euclidean Einstein equation whose boundary 
is the given Riemann surface; since none of the other fields are sourced by the boundary 
conditions, the solution will be locally AdS3 regardless of the matter content. In fact, we can 
go further: because the solution is locally AdS3, it will not be changed even in the presence 
of higher-derivative corrections to the bulk action, and the partition function will be changed 
only by an overall factor which amounts to a renormalization of the central charge. Such 
corrections correspond to moving away from infinite coupling in the boundary theory. Going 
even further, there is evidence that "free" large-c CFTs such as symmetric-product orbifolds 
also have the same partition functions as holographic ones [32]. (By "large-c theory" we 
mean one where the spectrum does not decompactify in the large-c limit, i.e. the number of 
operators below any given scaling dimension remains finite.) 

If all large-c CFTs have the same partition functions on arbitrary Riemann surfaces, then 
they also have the same Renyi entropies for arbitrary ./V and n. Further evidence for this 
proposition comes from several directions. First, the Ryu-Takayanagi formula [33, 34, 35] 
gives the same results for the entanglement von Neumann entropies of any set of intervals in 
any theory whose ground state is represented by AdS3 (global or Poincare), irrespective of 
what matter content the bulk theory might have. While this formula only applies when the 
bulk theory is Einstein gravity, it can be argued that this agreement survives higher-derivative 
corrections to the bulk action, based on the symmetries of AdS3 together with a standard 
ansatz for the effect of such corrections on the entropy [15]. Finally, in [15], direct evidence 
was found using CFT techniques that the Renyi entropies are the same for all large-c CFTs. 

We have argued that a large class of theories have identical partition functions and 
entanglement spectra. 12 All of these arguments are approximate, in that they apply only to 
the leading (order-c) parts of the partition functions and entropies. It seems likely that 1/c 
corrections will indeed distinguish between theories. 
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A. A resummation of the self-dual boson partition function 

In this appendix we will derive (3.23) from (3.20). Let = A + iK, with A,K symmetric 
and real, and K positive definite. Plugging (3.21) into the sum in (3.20), 



re 



(0|2fi)| 2 = X X exp [2m(n + efA^ 



-2vri(m + efA^m + e) j - 2vr(n + e) i K ij (n + e) j - 2vr(m + efK^im + e) j ] 

(A.I) 

Using 

2x l K ljX i + 2y^,y' = (x- yYKij (x - y) j + (x + yfK^ (x + y y , (A.2) 
the sum in (A.l) can be rewritten as: 



V 



1= ^2 ^ exp \2m(n + e) i A i j(n + e) j - 2iri(m + efA^m + e) 

-7r(n + m + 2e) i K ij (n + m + 2e) j - -rr(n - mfKij{n - m) j ] (A. 3) 

Defining £ = (n — m), we find n + fh + 2e = £ + 2m + 2e. Eq. (A. 3) can now be rewritten as: 

i= £ ex p 

-ir(£ + 2m + 2e) i K ij (l + 2m + 2e) 3 - irf*K ij e s ] (A A) 



-iri(2£ + 2m + 2e) i A ij (2£ + 2m + e) j - -m(2m + 2e) i A ij (2m + 2e) j 



m appears only in the form 2m + 2e. Summing over all m means summing over even 2m; 
since e G {0, summing over m and e means 2m + 2e take all values in Z ff once; we can 
thus replace the sums over rh and e with a sum over k = 2m + 2e, so that (A. 4) becomes 

I = Y exp [iri{2£ + k) i A ij {2£ + k) j - mtfAijW - ir(£ + kfKij^ + k) j - vrf K^ff] (A.5) 
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Shifting k — )■ k — £, (A. 5) becomes 

I = ^ exp •_>77/V',\ // /.' - TrfKijf - nk i K ij k j ] , (A.6) 

and recalling that the self-dual momentum lattice is simply ki±kji G \/2Z, we recover (3.23). 
B. Representations of the dihedral group 

The dihedral group D n is generated by two elements: g such that g n = 1, generating a Z„ 
subgroup, and r such that r 2 = 1, which satisfy the relation gr = rg^ 1 . A good reference for 
this subject is [36]. 

For n odd there are irreducible representations: 

• The trivial representation 1. 

• An additional one-dimensional representation T with g = 1, r = — 1. 
. n-l 



^— two-dimensional representations -R fc=1 r^i with 



\ (01 

e=^ 5 T " 10 



(B.l) 



One may write these representations for any k, but for k > (n — 1) /2, they are equivalent 
to one of the above. 

For n even there are irreducible representations: 

• The trivial representation 1. 

• The one-dimensional representation T with g = l,r = — 1. 

• The one-dimensional representation S with g = — l,r = 1. 

• The one-dimensional representation U with g = — l,r = —1. 

• ^^tj— two-dimensional representations i?. _ with 



e-^r \ (01 
e^ 5 5 T " 10 



(B.2) 



One may write these representations for any k, but for k > (n — 2) /2, they are equivalent 
to one of the above. 
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C. The boson at arbitrary radius 



Using the replica trick and classic results on correlators of twist fields in orbifold CFTs, 
Calabrese, Cardy, and Tonni (CCT) [6] calculated the Renyi entropies S^(A) for two intervals 
for the compact boson at arbitrary radius R; their result is shown in (2.4). In this appendix 
we will use the techniques developed in Section 3 to reproduce their result in a different way. 

The basic idea is do a direct calculation of the partition function on the replicated surface 
^n,2 by Weyl-transforming it to a non-singular surface and then applying (3.20) (or more 
precisely its generalization to arbitrary radius). However, to avoid having to compute the 
Liouville action (3.2) and the Laplacian determinant appearing in (3.20), we consider a ratio 
of the partition function of the boson to that of the Dirac fermion, and use Casini, Fosco, 
and Huerta's result (2.3) for the latter [8, 9]. All of the complicating factors are the same for 
the two theories, since they have the same central charge and same oscillators, and therefore 
cancel in the ratio, leading to quite a simple calculation. Furthermore, this derivation explains 
why the boson result takes the form of the fermion result plus a correction term. 

Let us proceed with the calculation. Since both theories have c = 1, they have the same 
Liouville action for any Weyl transformation; therefore the ratio of partition functions is 
independent of the choice of fiducial metric: 

s«(A) = 5® (A) + ^ ln W) = S - ){A) + rhi ln W) ■ (ai) 

The partition function for the boson at arbitrary R is essentially given by the same formula as 
in the self-dual case (3.23), but with the momentum lattice r( b ) (defined in (3.5)) generalized 
to 

r R = {(fc Lj k R ) : k L + k R e 72M, k L -k R e y^z} (c.2) 

(recall that rj = R 2 /R^ d ). It's useful to express k^^ in terms of the integral momentum and 
winding numbers n, w: 

kr p = n±\ -w. (C.3) 

The dot product appearing in the partition function becomes 

k ■ k' = k^k'r + k R k' R = —Tin' + rjww' . (C4) 

V 

As we showed in subsection 3.4, for N = 2, the period matrix is always imaginary (A = 0, 
$7 = iK). Hence the ratio of partition functions is 

z R _ 9 (i- fl )/2 tf(0|"/q)fl(0to") rr^ 

z(o " mm 2 ' 

As discussed at the end of subsection 3.3, the factor of 2^ 1-9 ^ 2 contributes a constant to the 
non-universal part of the Renyi entropies, and can be neglected. We thus have 

^)^) + ^.n «g« (C) 
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This result is extremely similar to the CCT result (2.4), the only difference being that the 
matrix T appearing there is replaced by the period matrix 0. However, CCT claim based on 
numerical evidence that for any rj, 

0(o|»/n) = tf(o|7?r) (c.7) 

(see appendices A, B, C of [6]), establishing that (C.6) and (2.4) agree. In fact, CCT give an 
explicit formula for Q, which is no more complicated than that for T, so (C.6) may itself be 
useful for direct calculation of the Renyis. 

Note that the result (C.6) holds not just for the case N = 2, arbitrary n, treated by 
CCT, but also for the case n = 2, arbitrary N, since f2 is imaginary there as well (as shown 
in subsection 3.4). However, to apply (C.6), one would need to compute £1. The Riemann 
surface T,2,n is hyperelliptic in that case, and one might be able to compute its period matrix 
via the associated Picard-Fuchs equations. 
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